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We analytically study gravitational radiation from corotating binary neutron stars composed of 
incompressible, homogeneous fluid in circular orbits. The energy and the angular momentum loss 
rates are derived up to the first post-Newtonian (1PN) order beyond the quadrupole approximation 
including effects of the finite size of each star of binary. It is found that the leading term of finite 
, size effects in the 1PN order is only 0(GM i ,/c 2 a i ,) smaller than that in the Newtonian order, where 

■ GM»/(?a* means the ratio of the gravitational radius to the mean radius of each star of binary, and 

the 1PN term acts to decrease the Newtonian finite size effect in gravitational radiation. 
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Around the beginning of the next century, the laser interferometers for detection of gravitational waves, such as 
LIGO §, VIRGO §, GEO600 §] and TAMA300 H) will be in operation. The most important targets for these 
detectors are coalescing binary neutron stars (BNS's) because in a period of their inspiraling phase they will emit 
gravitational waves of frequencies in the sensitive range of these detectors, i.e., from 10Hz to 1000Hz. We have to 
prepare accurate theoretical templates in order to extract informations of BNS's such as each mass, spin and so on 
. . | using matched filtering technique § . 

In their early inspiraling phase, the hydrodynamical effect of each star of binary is less effective and the point 
particle approximation works well. In the point particle approximation, the energy loss rate is calculated up to 2.5 
post-Newtonian (PN) order for arbitrary mass binaries |(|. The other approach for the inspiraling phase is the black 
hole perturbation method. In this formalism, the energy loss rate is calculated up to 4PN order in a Kerr black hole 
case H, and 5.5PN order in a Schwarzschild black hole case Isj. 

In the late inspiraling phase, however, the hydrodynamical effects of each star of binary become important, and we 
have to study the evolution of the binary taking into account them H jl(| |ll]] [|l2| . The importance of the study 
of the late inspiraling phase comes from the existence of the innermost stable circular orbit (ISCO). If we are able 
to know the location of the ISCO from the signal of gravitational waves, we can see a strong general relativistic 
phenomenon as well as we may get the information of the equation of state of neutron stars, i.e., the relation between 
the mass and the radius jl3) fllif . To know the precise location of the ISCO, it is desirable to construct the theoretical 
template of gravitational waves in the late inspiraling phase as in the early one. Up to the Newtonian order, we have 
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already understood gravitational radiation from binary systems of finite size stars |15 Jl6| |17| , but apparently the 
Newtonian treatment is not appropriate because BNS's are general relativistic objects. 

To investigate general relativistic effects on gravitational waves in the late inspiraling phase, in this paper, we 
analytically calculate the energy and angular momentum loss rates of the 1PN order including the finite size effects. 
In the calculation, we assume that BNS's are composed of the incompressible and homogeneous fluid and then obtain 
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an ellipsoidal equilibrium configuration for each star. Also, we suppose the corotating circular orbit in which each 
star of binary uniformly rotates around the center of mass of the system. Although the assumption of corotation is 
not appropriate for the realistic one which will be almost irrotational [Q p7j , the energy and angular momentum 
loss rates we calculate in this paper will be useful to see contribution of the finite size effects. 

This paper is organized as follows. In Sec. II, we derive the mass quadrupole moments and their time derivatives 
which are necessary for calculation of gravitational waves by using Blanchet-Damour formalism p9[ . In Sec. Ill, 
we give the rates of the energy and angular momentum loss gathering the time derivatives of the mass quadrupole 
moments which we derive in Sec. II, and their quantitative results are given in Sec. IV. Section V is devoted to 
summary and discussion. 

Throughout this paper, we use the unit G = 1 and c denotes the speed of light. Latin indices i, j, fc, • ■ • take 1 to 
3, and <5y denotes the Kronecker's delta. We use 7y and -fy as the quadrupole moment and its trace free part, 
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where M and denote the Newtonian mass and the principal axes of the star. 



II. FORMULATION 

In this paper, we calculate the energy and angular momentum loss rates of gravitational radiation from corotating 
binary stars. Each star of binary is assumed to be composed of the incompressible and homogeneous fluid because 
all the calculations can be done analytically. We pay particular attention to BNS's of equal masses in circular orbits 
whose coordinate separation is R. We assume that the center of mass of the binary system locates at the origin 
of the inertial coordinate system (Xi, X2, X3), and the center of mass of a star (star 1) locates at the origin of 
the comoving coordinate system {x\, X2, X3). In the comoving coordinate system, the other star (star 2) locates at 
(xi,X2,xs) = (— R, 0,0). The transformation from the coordinate system (Xi,X 2 ,X 3 ) to (xi, x 2 , £3) is given by 



X\ = ^xi + — ) < os ill .r_> sill < >/. 

X 2 = ^Xi + ^ ) sill*)/ ■+ .r-2 cos ill. 
X 3 = X 3 . 
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According to Blanchct and Damour |l9| , in the 1PN order, the energy and angular momentum loss rates are written 
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where My... and Sy are Blanchet-Damour's mass multipole moments and the current quadrupole moment, and the 
superscript (n) where n = 2, 3, 4 denotes the time derivative d n /dt n . Since we assume that the binary rotates around 
the X 3 axis, there is only i = 3 component in the angular momentum loss rate. In the subsequent subsections, we 
derive My... and S%j which are necessary for calculation of gravitational waves. 



A. Mass quadrupole moment 



First, we calculate My up to the 1PN order. My is the even term of the inertial coordinates Xi, X2, and A3, 
and we consider the equal mass binary. Then, the contribution from star 1 is the same as that from star 2. My is 
expressed as 
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In the following calculation, we integrate only over star 1 because the total mass quadrupole moment is twice as large 
as that of star 1. Moreover, as shown in Ref. p0| , the deformation from the Newtonian ellipsoid becomes higher order 
effect of R . Therefore, we integrate over the ellipsoid simply. 



1. Newtonian order 



The quadrupole moment of the Newtonian order appears in the first term of Eq. (2.4) and expressed as 

Dij = J d 3 Xp(x i X j - (2.10) 

Using Eq. fl2.l| ), we get the components as 

D n = D cos 2 tit + const., (2.11) 
D22 = -Dsin 2 tit + const., (2.12) 
D 12 = D sin tit cos tit, (2.13) 

where D = — I22 + R 2 M/A. In these equations, we use M = J pd 3 x as the Newtonian mass. Since D33 does not 
depend on time and the other terms vanish, we do not write them here. Note that in these components, only the 
contribution from star 1 is included. 

2. 1PN order 



Next, we calculate the 1PN order terms. For the 1PN order terms in Eq. (2.4), we define 

Vij = 2 f d 3 Xpv 2 X ih (2.14) 



Ulf 1 = 2 J cPXpU^Xij, (2.15) 
U 2 ^ 1 =2 f d 3 XpU 2 ^ 1 X ijl (2.16) 



Pij =3 / d 3 Xp-X l] , (2.17) 



Qy r = -— / d 3 XpX ij X 2 , (2.18) 
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because it is convenient to calculate the terms separately. In these equations, we use 
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where U 1 ^ 1 , U 2 ^ 1 and P denote the Newtonian potential generated by star 1 itself, the Newtonian potential generated 
by star 2, and the pressure, respectively fl(i|j . Ay... are index symbols introduced by Chandrasekhar [^(| and A = 
^2iAia 2 is calculated from EG] 



A = aia 2 a 3 



aia 2 a 3 



du 



o \J (a? + u){a\ + u)(a 2 + u) ' 
o ^(1 + i)(a 2 +i)(a§ +f) 



(2.24) 
(2.25) 



where 0:2 = 02/01 and «3 = 03/01 are axial ratios. 

When we derive the 1PN order terms, we need to calculate only some combinations: The 1PN order terms in 
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where we use the relation D 
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Then, we need to calculate only two combinatio ns of the components 



expressed as [(1, 1) — (2, 2)] and [(1, 2)]. We separately show the terms of the 1PN order in Eq. ( 2.26 ) in the following. 
(1.1) V u - V22: 

2Q 2 J d 3 Xp(X 2 + Xl){Xl - Xl) 

20 2 t>cos2m, (2.27) 



V u - V x 



where 



- _ R A M 3R 2 

v — 1 2~ 1 1111 — ^2222, 



(2.28) 



and 



Hjkl 



M 

35 



d X px^ X j X X 1 - 

3of (i=j = k = l) 

hi 



(2.29) 



4 



(1.2) V 12 : 
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= fi 2 0sin2m. (2.44) 



B. Mass octupole moment and current quadrupole moment 

In this subsection, we discuss on the mass octupole moment Mijk and the current quadrupole one Sij. They are 
written as 



where 



Myfc(t) - / d*XpX ijk , (2.45) 
Sij(t) = J d*Xpe kl<i X j>k v u (2.46) 



A< i B J> - i(A t B, + BiA,-) - i%(A fc B fc ). (2.47) 
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Mijk and Sy are the odd terms of the coordinates X±, X2, and X3. Therefore, the contribution from star 2 has the 
opposite sign of that of star 1. Then, the total and S%j vanish as 

(Myfcjtot = 1 M tjk + 2 M ijk = 0, (2.48) 
(S ij )tot = 1 S ij + 2 S ij = 0, (2.49) 

where in the right hand side, superscripts 1 and 2 denote the contributions from star 1 and star 2, respectively. 
Therefore, they have no contributions to the energy and the angular momentum loss rates in the case of the identical 
star binary. 



III. THE ENERGY AND ANGULAR MOMENTUM LOSS RATES 



Our purpose is to know the finite size effects of the 1PN order in gravitational radiation. In the following, we 
calculate the energy and angular momentum loss rates up to 0(R~ 3 ) in the 1PN order beyond the Newtonian order 
of the quadrupole formula. 



A. The energy loss rate 



Taking into account the fact stated in subsection II B , the energy loss rate is expressed as 
(*4 
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The contribution from the Newtonian order is 

(2D$)(2D$) = 128fl 6 D 2 
The contributions from the 1PN order are written as follows: 



Collecting the results we derived above, we can get the energy loss rate as 
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To express Eq. (|j|), we have used M, but it does not conserve through the sequence of the binary. Moreover, we 
simply used R as the orbital separation, but in the 1PN case, the center of mass deviates from that in the Newtonian 
one and as a result, the orbital separation also deviates from R. We need a conserved mass and an appropriate 
definition of the center of mass. 

If we use the conserved proper mass 
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and define the center of mass by it as 
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P* = P 

the orbital separation should be replaced by 
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Using these expressions, the energy loss rate is rewritten as 
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and 01*02*03* = M*/(47rp/3) = constant in the 1PN order. The expression of Eq. ( 3.15 ) seems to depend on the 
internal structure of the star of binary even if we take the limits di/R — > 0. Thus, Af* is not appropriate in comparing 
the energy luminosity with that for the point particle case. This is simply because Af* is not the gravitational mass. 
Instead of Af*, we need an appropriate gravitational mass. 
If we use the parameterized post-Newtonian (PPN) mass ]2l| 
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with the center of mass of each star defined as 
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and the orbital separation as 
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the leading inte rnal structure dependent term in the 1PN order is renormalized. Substituting the PPN mass and 
i?ppN into Eq. (|3.9|), we rewrite the energy loss rate as 
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Then we recover the energy loss rate given for the point equal mass binary case || when we take the limits a-i/R — » 0. 
(However, Mp pn a ppears to be the imperfect definition of gravitational mass as shown below.) 

From Eq. ( 3.20 ), we find that the leading finite size effect of 1PN order is only of 0(pAo)(^ 0(M/a,i)) smaller 
than the Newtonian finite size term. The next order terms are of 0(M/R) smaller than the Newtonian term. Here, 
we note that the effect of the spin-orbit (SO) coupling appears in the same order as the latter terms for the case of 
the corotating binary case. Including the SO coupling term, the energy loss rate of the point particle binary of equal 
mass is written as |22|| 
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We can see from this equation that the latter finite size terms of the 1PN order in Eq. (3.20) is not explained only by 
the SO coupling term. In order to examine whether this fact is really true or not, we rewrite dE/dt using the angular 
velocity which is the invariant value of the coordinate. 

Using the angular velocity, the orbital separation is written as 
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Also if we use Af*, the orbital separation is written as 
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Substituting these expressions into the equation of the energy loss rate, we have 
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On the other hand, in the point particle case, the energy loss rate is written as 
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Comparing the last term of Eq. (3.25) with that of Eq. (3.27), we find that they do not agree. A part of th e 
reason seems due to the fact that the 1PN quadrupole dependent term really appears in the last terms of Eq. ( (3.25 ) . 
However, we cannot still explain the disagreement completely because there remains terms except for the SO coupling 
one even if we take a spherical limit &ippn = a2PPN = 03PPN- Then, what is the reason why the terms remain? 
We guess that this is due to the imperfection of the definition of gravitational mass. We adopt the PPN mass as 
gravitational mass of each star to write the equation of the energy loss rate, but the PPN mass can renormalizc 
only the self-gravity and spin into mass. In the context of this paper, we need a conserved gravitational mass which 
renormalizes the quadrupole term produced by the tidal force of the companion star. As far as we know, however, 
no one has proposed such a mass. Nevertheless, -Mppn is not expected to be different from "true gravitational mass" 
M true so much because 
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Hence, we regard AfppN as an approximate gravitational mass in the following. 



B. The angular momentum loss rate 



Next, we calculate the angular momentum loss rate. As shown by Ostriker and Gunn p3[ , the angular momentum 
loss rate has the relation 



dJ 3 _ 1 dE 

~dT ~ n~dt 



(3.29) 

for any stable, uniformly rotating, equilibrium configuration. Then, the angular momentum loss rate is written as 



dJ 3 
~dT 



4M, 



PPN 



5c 5 



(2M PPN ft) 



7/3 



(2A/p PN ^) 4 / 3 



5Af PPN 



2 (a 



1PPN — a 2PPN) + 2a 2 ppN ~ a 2PPN 



t 3PPN 



373 
""84" 



(2M PPN ft) 2/3 
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525A/^ PN ( 2A/ppN ^) 4/3 | 160 ( Q lPPN ~ a 2PPN) + 77 ( 2a lPPN - a 2PPN - a 3PPN) 

( 840Mgpf ( 2602a iPPN - 1697 a 2 ppN - 443a 2 PPN ) 



(3.30) 



or 



dJ 3 
dt 



4M. 
' 5c 5 



(2M*Q) 7 / 3 



5M2 



2 - a 2 
2* a 3* 



>M»*_ 373 )2/3 

3 84 v ; 175M 2 



2 "'"' lo * (2M^\ 100(4 ~ 4) + 49(2a?„ " 4 - 4) 



(2M»Q) 2 



(2532a?, - 17674 - 4434) 



840M 2 

We will also show the direct calculation in Appendix A. 



(3.31) 



IV. NUMERICAL RESULTS 

In this section, we calculate the energy loss rate along the equilibrium sequence of the binary. The method is as 
follows. 

(i) Using the equilibrium equations for a corotating binary, we construct the sequence of the binary. The equations 
are written as 



- — = -wpa x Ai + — a t + — a u 



Pa 2 a M 2 

-_ = - 7rpa2 A 2 -—a 2 - 

P0 2 A M 2 



O 2 

5i N 2 

— a 2 , 



(4.1) 
(4.2) 
(4.3) 



where Qn denotes the Newtonian angular velocity written by 

2M 18^ii 



n 2 - 

~ R? 



R 5 



(4.4) 



At this stage, we determine R — R/a\, a 2 — a, 2 ja\ and = aa/a±. 

(ii) Inserting the equilibrium figures into the equation of the orbital angular velocity of 1PN order 



ft = 



( 2M* 
( R 3 



1/2 



3 2 2 2 \ 1 / *pA * . 11M, 

1 H ?r(2aT a , — a,,, — at A -< 

lOR 2y u 2 * 3 * J c 2 \ 5 4R 

29npAo* 2 2 2 \ | M * / 1t; , 2 71 2 
"( 2 «i* - a 2 * - a 3 J + — ^r(154a 1# - 71a 2>t 



50i? 2 



40i? 3 



65a q 



(4.5) 



and the equation of the total energy of the binary [Hj], we determine the initial and final points. Here, we regard 
the point which satisfies the condition fl/ir = 10Hz as the initial one and the point which satisfies fl/ir = 1000Hz or 
reaches the energy minimum (ISCO)| as the final one. 

(iii) Substituting the equilibrium figures into Eq. ( 3.25| ), we are able to have the energy loss rate. However, Eq. ( 3.25| ) 
is not written by conserved quantities. Then, we express Eq. (3.25) as 



Lb 



1 + i?N-fmite + E 



1PN 



E 



lPN-finitc 



(4.6) 



*In the corotating binary case, we call it the innermost stable corotating circular orbit (ISCCO) |10| . 
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where 



N 



dt I \ dt 
dE\ 2 



ic5 (2A/ PPN f!) 10 / 3 , (4.8) 



25n- 



N-finitc 



\ 1 

) (W7^{ 2(1 -^ ) + 2 -^-^}' (49) 



lPN-finitc 



1 /2M*ft\ 4/3 
5CfV^~J (aaos) 

2/3 

"84" V 

1 /2M* 0\ 4/3 A 



373 /2M. fi\ 

^ipn = -— (^2-7) > ( 4 - 10 ) 



350a 

In these equations, we use the compactness parameter 



C, S . (4-12) 



In calculation of the energy loss rate, we neglect the last term of Eq. ( 3.25| ), because this term includes the ambiguity 



of the definition of mass and we cannot separate the SO coupling term from it as mentioned in the last of subsection 



III A 



We repeat this procedure changing the compactness of the star M*/c 2 a* and the conserved proper mass M*. 
Figures l(a)-3(c) show the results we have using this procedure. It is found from these figures that the finite size 
effect of the 1PN order is only a few factor less than that of the Newtonian order. We can explain this feature by 



comparing Eq. (|4.9D with Eq. (4.11). Essentially, the difference between these equations is 0{C S ). Therefore, the 



finite size effect of the 1PN order becomes 0{C S ) less than that of the Newtonian order. 



There is another feature to mention specially. Comparing Eq. ( |4.9| ) with Eq. (4.11), it is found that the contribution 
from the 1PN order has the opposite sign of that from the Newtonian order. We explain this fact as follows: Including 
the 1PN order terms, the self-gravity of each star of binary becomes stronger than that in the Newtonian gravity. 
Then, the stars become more compact, and it is more difficult to deform them. This leads to the decrease of the finite 
size effect. Moreover, it is found from these figures that when we increase the compactness parameter fixing M*, both 
the finite size effects of the Newtonian order and the 1PN one decrease, and also the difference between their absolute 
values decreases. The explanation for these behaviors is the same as the above one, i.e., the self-gravity of each star 
becomes stronger. 

Finally, we discuss on the inclinations of the lines in the figures. First, the inclination of the 1PN line becomes 2/3 
in the log- log plot figures because of Eq. ( 4.10| ). Next, we can see from the figures that the inclination of the finite 



size effect lines is ~ 3.4. This is because the deviation from the spherically symmetric star is made by the tidal force 



M/R 3 . We find from Eq. (3.1C) that M/R 3 is equal to ft 2 . Then, combining this fact with Eq. flP|) or (4.11), we 



can conclude the inclination of the finite size effect lines becomes ~ 10/3. 

In Table I, we show the results of the initial and final orbits. Note that the frequencies of the final orbits are 1000Hz 
in the cases of (Af*/M Q , M„/c 2 a») = (1.4, 0.20), (1.4, 0.25), (1.6, 0.25), and (1.8, 0.25), while in the other cases, 
the ISCCOs are the final orbits. Also we present the PPN mass (gravitational mass) at infinity which is written as 

Mp PN , inf = M* M - 3 -C s ) (R -» 00). (4.13) 



V. SUMMARY AND DISCUSSION 



In this paper, we have analytically calculated the energy and the angular momentum loss rates. The conclusions 
are as follows. 

(i) The leading finite size term of the 1PN order in dE/dt reduces (ciEJ/e^N-fimte by ~ 40% which is 2C S times as 
large as that of the Newtonian order. 

(ii) The reason of the weak convergence of the finite size terms is that both of them are in the same order when we 
see Eq. Q3.25 ) in the viewpoint of the series of ft. However, the 1PN term is of 0(C S ) higher than the Newtonian one 
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in the viewpoint of the PN approximation. Therefore, the 1PN term is only by a factor of 0(C S ) smaller than the 
Newtonian one. 

(hi) The leading finite size term of the 1PN order has the opposite sign of that of the Newtonian order. This is because 
the 1PN order terms make the self-gravity of the star stronger, then it is difficult to deform the star. Therefore, the 
finite size effects on the energy loss rate decrease. 

Next, we discuss the definition of gravitational mass of each star. In this paper, we use the PPN mass as the 
gravitational mass in order to compare the equations of the energy loss rate with those of the point particle case. 
However, the PPN mass is appropriate only for the point particle case (i.e., as long as we ignore the tidal forces on 
each star of binary) pl| . In the context of present paper, we should use a gravitational mass which conserves and 
appropriately renormalizes the effect of the finite size of stars even if the tidal force exists. Unfortunately, no one has 
proposed such a mass as far as we know. To define such gravitational mass will be an unresolved problem. 

Finally, we point out a problem with regard to the back reaction by gravitational radiation as follows. In this paper, 
we have calculated the total gravitational wave luminosity. However, it is not trivial what fraction of it contributes to 
the back reaction to SI because the structure of each star (i.e., in the incompressible case) is also affected due to the 
radiation reaction. This situation is in contrast to that for the point particle case where we only need to consider the 
radiation reaction to f2. To know the change rate of f£ and ai in the 1PN order, we need to solve the hydrodynamic 
equation in the 1PN order including the 3.5PN radiation reaction terms. Here, we note that from the observational 
point of view, we need not the total luminosity, but the back reaction of i\. Hence, to resolve this problem seems one 
of the important issues in order to study the late time evolution of BNS's. 
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APPENDIX A: THE DIRECT CALCULATION OF THE ANGULAR MOMENTUM LOSS RATE 

In the case of the identical star binary, Myt and SL- are zero. Therefore, the angular momentum loss is written as 



d£_ 
~~dt 



2 

be 5 



(Al) 



For the quadrupole moment, the contributions from star 2 are the same as those from star 1. Then, the total 
quadrupole moments double those of star 1. 

In the case we take in this paper, the binary rotates around the X% axis. Then, the angular momentum loss has 
only i = 3 component and is written as 



dt 



5c 5 



{2Df)(2D 



1:1 



+ ^U 2 ^ ) )( 2 n 



,(3) 



1^1(3) 



2U, 



2-a(3) 

kl 



2P, 



(3) 



2Q 



(3) 
kl 



2R 



+(2D$) ( 2vf + 2V)r m + 2U"r 1(2) + 2/f } + 2Q« + 2R?> 



When we use the explicit forms for the quadrupole moments, the angular momentum loss is rewritten as 



dJ 3 
~dT 



8 

5? 



9 n(2)n(3) ,n{2) n (3) 



+(Pu - P 22 ) C3) + (Qn - Q 22 ) (3) + (Ru - i? 22 ) (3) ) 



(A2) 
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+ p(3) +Q (3) +i2 (3) 



l 12 



+(Pn - J p 22 ) (2) + (Qn - g 22 ) (2) + (fin - i? 22 ) (2) ) 

-2D<»> f V« + C/ 1 V 1(2) + ^V 1(2) + P[f + Q® + 



(A3) 



where we use the relation D 22 = -i) n '. 

The contribution from the Newtonian order is 



(A4) 



The contributions from the 1PN order are written as follows: 



-D$(Vu-V 22 )W+2D?M2 ) 



D$( Vll -V 22 )W-2D<$vW 



32fTW, 



= 32fl 5 DU 1 ^ 1 , 



(2) rr 2-l(3) 
12 



-D$(P n -P 22 )W+2D$P$ 



DY 2 >(Qn-Q 22 )W + 2D { $Q 



(2) n (3) 

12 



-^(fi 11 -fi 2 2) (3) + 2i?f 1 ) J R (3) 
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D? 2 \ulr X -Ut 2 ^ -2D^ul 2 m 
32n 5 DU 2 ^\ 

48f7 5 L>P, 

^ 3 2 ) (Q 1 i- 22 ) (2) - 2^?Q(f 

--n 7 DQ, 

D$(R n -R 22 )W-2D$R$ 
21 ^ 



Gathering these terms, we have the angular momentum loss rate as 



(A5) 
(A6) 
(A7) 
(A8) 

(A9) 
(A10) 



di 5c 5 N 



i + ^2(«?-«i) + o(i?- 4 ) 



+- 



113M TtpAo J 296 . 2 2x l/ 2 2 2n\ 

+ 84W (358a? + 161a2 ~ 2465a ' } + ° ( ^ 4) 



4M /2M\ 



7/2 



1 + i^i 5 (a ' ~ ° 2) + 2 { 1 - a 2- a t)}+ 0(R- 



113M ttjoAc 



168.R i? 2 



296 2 2 133 2 2 2x\ 

— (a! - a 2 ) + -rj(2«i - oa - a 3 ) j 



M 



1680i? 3 



(2674a 2 - 2017a^ + 3235a 2 ;) + 0(iT 4 ) 



(All) 



(A12) 



As in the case of the energy loss rate, there remains the internal structure dependent term in the limits ai/R^ 
when we rewrite the angular momentum loss rate using M* as 
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dt 



4M* /2Af : 



5c 5 V R 



7/2 



0(R~ 4 ) 



9npA 1285M* tt P A f 1672 



ifli - a 2 ) + y ( 2a i - a 2 - «3) 



840i? 3 



84i? 5i? 2 [ 105 
(64274a? - 41171o| - 19645a 2 ,) + 0(iT 4 



l.U. /2A£ 



7/2 



5 



^2 ^-u* ' 2 

_ y 

84i?» ~ 5.R 2 | 105 



QirpA 0te 1285A/* 7rpA * [3016 2 2 53 , 

-{a u - a 2 J + y(2a ll( , - a 2 „ - a 3 J 



840i?2 



(74998a 2 „ - 46813a 2 ,* - 22375a 2 ,*) + 0(iC 4 ) 



However, when we substitute the PPN mass into Eq. (A12), we have 

1/8/2 - .3 



dJ 3 
dt 



4Af PPN {2M PPN \ 7/2 



5c 5 V # 



i? 2 [5 



a z 2 ) + -{2ai-4-ai)\+0(R- i ) 



1 



1285M PPN vrpA f 32 



84i? 

A/ppn 

168.R 3 



" ( A A\ I /o A A A 

— [a x - a 2 ) + [2a x -a 2 - a 3 



5R 2 21 



(13006a 2 - 8083a| - 3929a 2 .) + 0(ir 4 ) 



4Af PPN / 2 M] 



PPN 



5c 5 V -Rppn 



7/2 



p2 

n PPN 



/ 2 2 \ /o 2 2 2 \ 

-(a 1PPN — g^ppn) + 2(2a 1PPN — a 2PPN — a 3PPN ) 



+0(R P p N ) 



1285A/ppn tt/oAoppn f 256 . 2 
^tt; rr^ < — a, 



84i?ppN 

A/ppN 

840i? PPN 



5i?p PN 



I 21 



1PPN ~ "2PPNV 



11 (2a 



1PPN — a 2PPN a 3PPN 



(75754a 2 PPN - 46057a 2PPN - 22375a 2 p PN ) + 0(i? PPN ) 



and this expression does not depend on the internal structure of the star. 



(A13) 



(A14) 



(A15) 



(A16) 



If we subst itute Eqs. (3.23) and/or (3.24) into Eqs. (A15) and/or (A1S), we have the same equations as (3.3C) 
and/or (3.31). 
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TABLE CAPTIONS 



Table I. The coordinate orbital separation R*/a*, axial ratios 0,2/0,1 and 03/01, and the frequency of gravitational 
wave f2/7r[Hz] of the initial and final orbits. 





Initial 


Final 


M*/M 




Mpp N ,inf/M 


R* / a* 


a 2 /ai 


a 3 /ai 


q/ty[Hz] 


R*/a, 


a 2 /ai 


a 3 /ai 


fi/7r[Hz] 




0.15 


1.274 


48.9478 


0.999983 


0.999971 


10.00 


2.6511 


0.888010 


0.837826 


706.08 


1.4 


0.20 


1.232 


62.9341 


0.999993 


0.999989 


10.00 


2.6278 


0.918423 


0.877899 


1000.0 




0.25 


1.190 


75.2153 


0.999997 


0.999995 


10.00 


3.1321 


0.965573 


0.945395 


1000.0 




0.15 


1.456 


44.7585 


0.999977 


0.999962 


10.00 


2.6511 


0.888010 


0.837826 


617.82 


1.6 


0.20 


1.408 


57.5478 


0.999991 


0.999986 


10.00 


2.4554 


0.898344 


0.851178 


959.53 




0.25 


1.360 


68.7777 


0.999996 


0.999993 


10.00 


2.8328 


0.953794 


0.927832 


1000.0 




0.15 


1.638 


41.3604 


0.999971 


0.999952 


10.00 


2.6511 


0.888010 


0.837826 


549.17 


1.8 


0.20 


1.584 


53.1789 


0.999989 


0.999982 


10.00 


2.4554 


0.898344 


0.851178 


852.92 




0.25 


1.530 


63.5560 


0.999995 


0.999992 


10.00 


2.5904 


0.939745 


0.907524 


1000.0 



TABLE I. 
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FIGURE CAPTIONS 



Figs.l. The energy loss rates relative to that of the Newtonian order are shown as a function of the frequency of 
gravitational waves /gw — fi/7r. The sequences are calculated from the initial orbit (10Hz) to the final one (1000Hz 
or ISCCO). In these figures, we fix the conserved proper mass as M* = 1.4M and change the compactness parameter 
as C s = 0.15 (a), 0.2 (b), 0.25 (c). Solid, dotted and dashed lines denote contributions from the finite size effect of 
the Newtonian order, that of the 1PN order and the 1PN order terms, respectively. The long-dashed line denotes the 
contribution from the Newtonian order term which is unity in these figures. In fig. 1(a), we present the location of the 
energy minimum point (ISCCO) using the solid vertical line. 

Figs. 2. The conventions are the same as those in figs.l, but for M*. Here, we fix it as M* = I.QMq and change the 
compactness parameter as C s — 0.15 (a), 0.2 (b), 0.25 (c). In fig. 2(a) and 2(b), we present the location of the energy 
minimum point (ISCCO) using the solid vertical line. 

Figs. 3. The conventions are the same as those in figs.l, but for M*. Here, we fix it as = 1.8M Q and change the 
compactness parameter as C s — 0.15 (a), 0.2 (b), 0.25 (c). In fig. 3(a) and 3(b), we present the location of the energy 
minimum point (ISCCO) using the solid vertical line. 
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Fig. 1(a) 
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Fig.l(b) 
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Fig. 1(c) 
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Fig.2(a) 
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Fig.2(b) 
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Fig.2(c) 



23 




Fig.3(a) 
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Fig.3(b) 
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Fig.3(c) 
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